Hydrodynamics of Suspensions of Passive and Active Rigid Particles: A Rigid Multiblob Approach

We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally-resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. For unbounded suspensions and suspensions sedimented against a single no-slip boundary, we rely on existing analytical expressions for the Rotne-Prager tensor combined with a fast multipole method or a direct summation on a Graphical Processing Unit to obtain an simple yet efficient and scalable implementation. For fully confined domains, such as periodic suspensions or suspensions confined in slit and square channels, we extend a recently-developed rigid-body immersed boundary method to suspensions of freely-moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid body equations converges in a bounded number of iterations regardless of the system size. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall.Comment: Under revision in CAMCOS, Nov 201

A late medieval warm period in the Southern Ocean as a delayed response to external forcing?

International audienceOn the basis of long simulations performed with a three‐dimensional climate model, we propose an interhemispheric climate lag mechanism, involving the long‐term memory of deepwater masses. Warm anomalies, formed in the North Atlantic when warm conditions prevail at surface, are transported by the deep ocean circulation towards the Southern Ocean. There, the heat is released because of large scale upwelling, maintaining warm conditions and inducing a lagged response of about 150 years compared to the Northern Hemisphere. Model results and observations covering the first half of the second millenium suggest a delay between the temperature evolution in the Northern Hemisphere and in the Southern Ocean. The mechanism described here provides a reasonable hypothesis to explain such an interhemipsheric lag

Mutations in Hydin impair ciliary motility in mice

Chlamydomonas reinhardtii hydin is a central pair protein required for flagellar motility, and mice with Hydin defects develop lethal hydrocephalus. To determine if defects in Hydin cause hydrocephalus through a mechanism involving cilia, we compared the morphology, ultrastructure, and activity of cilia in wild-type and hydin mutant mice strains. The length and density of cilia in the brains of mutant animals is normal. The ciliary axoneme is normal with respect to the 9 + 2 microtubules, dynein arms, and radial spokes but one of the two central microtubules lacks a specific projection. The hydin mutant cilia are unable to bend normally, ciliary beat frequency is reduced, and the cilia tend to stall. As a result, these cilia are incapable of generating fluid flow. Similar defects are observed for cilia in trachea. We conclude that hydrocephalus in hydin mutants is caused by a central pair defect impairing ciliary motility and fluid transport in the brain

A simple catch: Fluctuations enable hydrodynamic trapping of microrollers by obstacles

It is known that obstacles can hydrodynamically trap bacteria and synthetic microswimmers in orbits, where the trapping time heavily depends on the swimmer flow field and noise is needed to escape the trap. Here, we use experiments and simulations to investigate the trapping of microrollers by obstacles. Microrollers are rotating particles close to a bottom surface, which have a prescribed propulsion direction imposed by an external rotating magnetic field. The flow field that drives their motion is quite different from previously studied swimmers. We found that the trapping time can be controlled by modifying the obstacle size or the colloid-obstacle repulsive potential. We detail the mechanisms of the trapping and find two remarkable features: The microroller is confined in the wake of the obstacle, and it can only enter the trap with Brownian motion. While noise is usually needed to escape traps in dynamical systems, here, we show that it is the only means to reach the hydrodynamic attractor.“la Caixa” Foundation (ID 100010434), fellowship LCF/BQ/-PI20/11760014. The European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement no. 847648

Diffusive limits on the Penrose tiling

In this paper random walks on the Penrose lattice are investigated. Heat kernel estimates and the invariance principle are shown

High density plasma deposition of device quality silicon nitride. II. Effects of thickness on the electrical properties

Dielectric behavior of SiNx films, fabricated by microwave electron cyclotron resonance discharge, has been studied as a function of film thickness on the basis of the current–voltage and the capacitance–voltage characteristics. In the thickness range (20 nm,d,80 nm), the resistivity and the critical field for SiNx were found not to be sensitive to the film thickness ~d! and which was opposite to strong dependence of the dynamic dielectric constant e d on thickness. To explain the e d behavior as a function of d, a model based on trapped space charge effects is proposed. The dominant mode of electronic conduction, determined from J –E1/2 curves and Arrhenius plots of leakage current, appears to be Poole–Frenkel emission only for thicker films (d.20 nm). Finally, the spatial profile of fixed charges reveals that SiNx /Si interface has a much greater concentration of defects than the bulk film

What Are the Ants Doing? Vision-Based Tracking and Reconstruction of Control Programs

©2005 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other users, including reprinting/ republishing this material for advertising or promotional purposes, creating new collective works for resale or redistribution to servers or lists, or reuse of any copyrighted components of this work in other works.Presented at the 2005 IEEE International Conference on Robotics and Automation (ICRA), 18-22 April 2005, Barcelona, Spain.DOI: 10.1109/ROBOT.2005.1570762In this paper, we study the problem of going from a real-world, multi-agent system to the generation of control programs in an automatic fashion. In particular, a computer vision system is presented, capable of simultaneously tracking multiple agents, such as social insects. Moreover, the data obtained from this system is fed into a mode-reconstruction module that generates low-complexity control programs, i.e. strings of symbolic descriptions of control-interrupt pairs, consistent with the empirical data. The result is a mechanism for going from the real system to an executable implementation that can be used for controlling multiple mobile robots

Climate dependent contrast in surface mass balance in East Antarctica over the past 216 ka

ABSTRACTDocumenting past changes in the East Antarctic surface mass balance is important to improve ice core chronologies and to constrain the ice-sheet contribution to global mean sea-level change. Here we reconstruct past changes in the ratio of surface mass balance (SMB ratio) between the EPICA Dome C (EDC) and Dome Fuji (DF) East Antarctica ice core sites, based on a precise volcanic synchronization of the two ice cores and on corrections for the vertical thinning of layers. During the past 216 000 a, this SMB ratio, denoted SMBEDC/SMBDF, varied between 0.7 and 1.1, being small during cold periods and large during warm periods. Our results therefore reveal larger amplitudes of changes in SMB at EDC compared with DF, consistent with previous results showing larger amplitudes of changes in water stable isotopes and estimated surface temperature at EDC compared with DF. Within the last glacial inception (Marine Isotope Stages, MIS-5c and MIS-5d), the SMB ratio deviates by up to 0.2 from what is expected based on differences in water stable isotope records. Moreover, the SMB ratio is constant throughout the late parts of the current and last interglacial periods, despite contrasting isotopic trends.This is the accepted manuscript. It is currently embargoed pending publication

The influence of the synoptic regime on stable water isotopes in precipitation at Dome C, East Antarctica

Abstract. The correct derivation of paleotemperatures from ice cores requires exact knowledge of all processes involved before and after the deposition of snow and the subsequent formation of ice. At the Antarctic deep ice core drilling site Dome C, a unique data set of daily precipitation amount, type, and stable water isotope ratios is available that enables us to study in detail atmospheric processes that influence the stable water isotope ratio of precipitation. Meteorological data from both automatic weather station and a mesoscale atmospheric model were used to investigate how different atmospheric flow patterns determine the precipitation parameters. A classification of synoptic situations that cause precipitation at Dome C was established and, together with back-trajectory calculations, was utilized to estimate moisture source areas. With the resulting source area conditions (wind speed, sea surface temperature, and relative humidity) as input, the precipitation stable isotopic composition was modeled using the so-called Mixed Cloud Isotope Model (MCIM). The model generally underestimates the depletion of 18O in precipitation, which was not improved by using condensation temperature rather than inversion temperature. Contrary to the assumption widely used in ice core studies, a more northern moisture source does not necessarily mean stronger isotopic fractionation. This is due to the fact that snowfall events at Dome C are often associated with warm air advection due to amplification of planetary waves, which considerably increases the site temperature and thus reduces the temperature difference between source area and deposition site. In addition, no correlation was found between relative humidity at the moisture source and the deuterium excess in precipitation. The significant difference in the isotopic signal of hoarfrost and diamond dust was shown to disappear after removal of seasonality. This study confirms the results of an earlier study carried out at Dome Fuji with a shorter data set using the same methods